Published in: M. Nakagawa and S. Luding (eds.), Powders and Grains 2009: Proceedings of the 6th International 
Conference on Micromechanics of Granular Media, American Institute of Physics (Melville, New York, 2009), p. 859-862 



Fractal Substructures due to Fragmentation and 

Reagglomeration 

D. E. Wolf, T. Poschel^ T. Schwager**, A. Weuster* and L. Brendel* 

* Department of Physics and CeNIDE, University of Duisburg-Essen, 47048 Duisburg, Germany 
^CBI and Center of Excellence EAM, Universitdt Erlangen-NUrnberg, 91058 Erlangen, Germany 
**Charite, Augustenburger Platz 1, 13353 Berlin, Germany 

Abstract. Cohesive powders form agglomerates that can be very porous. Hence they arc also very fragile. Consider a process 
of complete fragmentation on a characteristic length scale i, where the fragments are subsequently allowed to settle under 
gravity. If this fragmentation-reagglomeration cycle is repeated sufficiently often, the powder develops a fractal substructure 
with robust statistical properties. The structural evolution is discussed for two different models: The first one is an off-lattice 
model, in which a fragment does not stick to the surface of other fragments that have already settled, but roUs down until 
it finds a locally stable position. The second one is a simpler lattice model, in which a fragment sticks at first contact with 
the agglomerate of fragments that have already settled. Results for the fragment size distribution are shown as well. One can 
distinguish scale invariant dust and fragments of a characteristic size. Their role in the process of structure formation will be 
addressed. 
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INTRODUCTION 

Most particles attract each other, be it by van-der-Waals 
forces or by microscopic liquid menisci at the con- 
tacts, but often this attraction is so weak compared to 
other forces acting on the particles that it may safely 
be neglected, when explaining the physical behaviour of 
granular matter. There are notable exceptions, however, 
where attractive forces are decisive. In this paper we con- 
sider one such case, in which the attractive force between 
two particles is much larger than their weight. This is true 
for nanopowders [1], but also for larger particles under 
microgravity [2], or for wet sand [3]. 

Such a granular medium with adhesion between the 
particles can sustain very high porosity in spite of grav- 
ity. This can be easily shown: Take two equal glass con- 
tainers, one empty and the other one at most half filled 
with dry sand. Porosity is known to be about 0.36 (ran- 
dom dense packing). Then about 10 volume percent of 
water is kneaded into the sand. One gets a smooth dough. 
By and by this dough is fragmented by means of a fork 
into crumbs which are poured into the empty container. 
One finds the filling height increased. Then the crumb as- 
sembly is again fragmented with the fork and filled back 
into the original container. Again the filling height in- 
creases. Repeating this procedure a few times leads to 
a sand packing with porosity of about 70 percent. Ther- 
mophoretic aerosol deposits can even have porosities as 
high as 99 percent. 

In this paper we report on the structural properties 
of a porous assembly of adhesive particles obtained by 



repeated fragmentation and reagglomeration. We show 
that a fractal substructure forms and that the fragment 
size distribution is very broad. To obtain statistically 
significant results for the structure one has to con- 
sider systems with more than one million particles for 
many fragmentation-reagglomeration cycles. This is be- 
yond the capability of Molecular Dynamics simulations. 
Therefore, two simplified two-dimensional models were 
studied. 

The first one, an off-lattice model [4], is a general- 
ization of a model for the sequential deposition of non- 
adhesive spherical particles under the influence of grav- 
ity [5, 6, 7, 8, 9, 10]. The second one is a lattice model, 
generaUzing the model of ballistic deposition [II, 12]. In 
both models the following procedure is repeated many 
times: First the agglomerate is cut with a square mesh 
into portions. The linear mesh size ^ can be viewed as the 
typical scale of the fragmentation process. A portion may 
consist of several disconnected fragments. The models 
differ in the way these fragments then settle under grav- 
ity. They do so as rigid bodies without taking adhesion 
forces with other particles into account. This is justified, 
if the flakes are sufficiently large, so that their weight 
exceeds the adhesive force between the particles. After 
this reassembly of the fragments the agglomerate is cut 
again with the square mesh, and so on, see Fig. I and 
Fig. 6. In the following, lengths are given in units of the 
average particle radius (off-lattice model) respectively 
the lattice constant (lattice model), masses in units of 
the particle mass, and time as number of fragmentation- 
reagglomeration cycles. 
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FIGURE 2. The asymptotic filling heigth Aoo grows as a 
power law h„,{£) ~ i" with mesh size £. The full line shows the 
best fit, a = 0.327. Inset: The relaxation time no{E) increases 
linearly with mesh size. 



FIGURE 1. Evolution of the packing in the off-lattice model, 
a) Initial packing generated by random sequential sedimenta- 
tion [5]. The packing is cut by a square mesh into fragments 
(£ = 20); b) The fragments are considered as rigid bodies and 
deposited (1" generation). Again the packing is cut by the 
square mesh (here not shown); c) the fragments are deposited 
again (2"'' generation), and so on; d) 3"^ generation; e) 4''' 
generation; f) 120'* generation. 



THE OFF-LATTICE MODEL [4] 

In this model the reagglomeration of the fragments is 
simulated in the following way: Each fragment starts 
at a random position with a random orientation well 
above the already deposited material (the configuration 
of which is regarded as frozen in). Following gravity, 
it moves downwards until it touches the bottom of the 
container or contacts another already deposited particle. 
Then it rolls down as a rigid body, again following grav- 
ity, until the vertical projection of its center of mass falls 
in between two points of contact (with the container or 
previously deposited particles). 

We simulated up to 3 million particles represented 
by discs with a narrow size distribution (10% variance). 
The initial state was a densely packed agglomerate. 
In each fragmentation-reagglomeration cycle the filling 
height increases. Asymptotically, the powder adopts a 
very porous, statistically invariant structure, which is ro- 
bust with respect to fragmentation at scale £ and subse- 
quent gravitational settling of the fragments. 

The filling height /!„ at iteration step n approaches the 
asymptotic height, exponentially with a relaxation 
time, HQ. The inset of Fig. 2 shows that 
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FIGURE 3. Average fragment mass as a function of the mesh 
size £ in the steady state of the off-lattice model. 



For the asymptotic filling height, a power law 

h^{£)ocf with a = 0.327 (2) 

gives a very good fit (see Fig. 2). This implies that the 
number of portions cut from the steady state configura- 
tion of a system of width W scales like A^p = h^W /f' °^ 
£"^^. Consequently, the mass of a £ x ^-portion has a 
fractal dimension d{, 



M 



■■ f 'f with df^2 



a 



1.67±0.03. (3) 



In the steady state the number of fragments per portion 
is determined by the fact that the fragmentation process 
cuts on average as many particle contacts as will be 
reestablished by the agglomeration process. Since every 
fragment, when settling, creates two new contacts, the 
number of fragments per portion must be proportional to 
the number of contacts cut at the boundary of a mesh cell. 
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FIGURE 4. Normalized fragment mass distribution for dif- 
ferent mesii sizes i. f(m) is the number of fragments of mass 
m divided by the total number of fragments for a given £. Inset: 
Data collapse using °= with df = 1.695, and T = 1.41. 
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FIGURE 5. Chunk mass nic as a function of mesh size. Slope 
of straight Une is 1.695. 



This is confinned by Fig. 3, which shows that M/Nf oc £. 

A detailed understanding of the fragment properties is 
provided by the distribution of fragment masses, shown 
in Fig. 4. Two types of fragments must be distinguished, 
large chunks at the upper end of the mass spectrum 
with a characteristic size mc and scale invariant dust 
responsible for the power law part that is cut off by 
OTc- Comparing the mass distributions for different mesh 
sizes i shows, that they can approximately be written in 
the form 

^('"'^)='""^^~(^)' 

where the scaUng function f{x) is constant for x <i; 1, 
goes through a maximum at x = 1, and has an approx- 
imately Gaussian tail for x 1. The typical mass nic 
of the chunks has a power-law dependence on the mesh 
size, rric = 0.304 f}-^'^^ (Fig. 5), the exponent being in 
good agreement with the value of df, Eq. (3). 

The fractal chunks are only a tiny fraction of the total 
number of fragments, A^f. This follows from the fact that 
the width of the chunk-distribution is proportional to rric. 
Therefore, the number of chunks, A^c. divided by the 
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FIGURE 6. Evolution of the packing in the lattice model. 
a)l'" generation. The initial packing was cut by the indicated 
square mesh into fragments {£ = 32); b) 4''' generation; c) 
10''' generation; d) 20'* generation; e) 40'^" generation; f) 100''' 
generation. 



number of fragments, Nf, decreases with increasing mesh 
size Uke 



»df(I-T) 



(6) 



The fractal dimension of the chunks and the dust ex- 
ponent T are not independent of each other but obey the 
scaling relation 

rff(2-T) = l. (7) 

The reason is the following: As the chunk mass mc °^ (f'^ 
scales in the same way as the total mass per portion, 
M/Np, the number of chunks per portion cannot depend 
on £. Using (4) and (6) one concludes that 



const. : 
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. £df(l-T)^(if- 
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which proves (7). For df = 1.695 this implies T = 1.41. 
These values lead to an excellent data collapse for the 
fragment mass distributions (see Fig. 4 (inset)). 

We have seen, that the overwhelming number of frag- 
ments are dust particles, apart from a vanishing fraction 
(6) of chunks. This explains, why the mass (essentially 
mass of chunks) per fragment (essentially per dust parti- 
cle) has nothing to do with the fractal dimension, but is 
proportional to £ (see Fig. 3) . 



THE LATTICE MODEL 

In order to explore how universal the fractal substructure 
is, we studied another fragmentation-reagglomeration 
model. Here all particles hve on a square lattice (of width 
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FIGURE 7. The steady state bulk density v as a function of 
the mesh size £ for the lattice model. 



W and height H) with lateral periodic boundary condi- 
tions, which in each cycle is subdivided into square cells 
of size £ X £. Then, traversing from bottom to top and 
each row from left to right, the content of each cell is con- 
sidered. Disconnected fragments therein are identified 
and - without changing their orientation - subsequently 
deposited into an (initially empty) W x //-container at a 
randomly chosen horizontal position. Deposition means 
here, that the fragment sinks vertically until it forms the 
first vertical contact with the container bottom or a pre- 
viously deposited cluster. At this point it stops moving. 
In contrast to the ballistic deposition model [12] horizon- 
tal contacts have no effect for the deposition process, but 
they become as sticky as any other contact, once the frag- 
ment is at rest. As the fragments do not roll down until 
they form a second contact, the structure, see Fig. 6, is 
more treelike than in Fig.l. 

In the following, we show results for a system with 
N = 2^^ particles (occupied lattice sites) in a container 
of width W = 4096. The mesh size ranges from £ — 
16, . . . ,256. Data are averaged over 6 independent runs, 
starting from initial conditions of randomly deposited 
single particles. As for the off-lattice model the system 
approaches a steady state in relaxation time «o °= ^, i-e. 
the exponent z is 1 . 

Figure 7 shows that the bulk volume fraction v in 
the steady state is a power law of the mesh size £ with 
a non-trivial exponent: v oc £d-f-'^ = £-0.387±0.005 n^^^ 
indicates again a fractal substructure up to scale £ which 
is confirmed in Fig. 8. The fractal dimension df^ 1.5 ± 
0.1 obtained from the box counting method should be 
more reliable than the one extracted from the global solid 
fraction. Fig. 7, because of the crossover from the fractal 
structure on small scales to the homogeneous density 
on scales larger than £. For the lattice model the fractal 
dimension seems to be a bit smaller than for the off- 
lattice model. Lattice anisotropy has a similar effect on 
the fractal behaviour of DLA [13]. 



FIGURE 8. The box counting method (e being the box size) 
applied to the steady state bulk structure of the lattice model: 
Below the mesh size I, a fractal dimension of about 1.5 prevails. 
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